Mean Field Theory for Pedestrian Outflow through an Exit 
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An average pedestrian flow through an exit is one of the most important index in evaluating pedes- 
trian dynamics. In order to study the flow in detail, the floor field model, which is a crowd model 
by using cellular automaton, is extended by taking into account a realistic behavior of pedestrians 
around the exit. The model is studied by both numerical simulations and cluster analysis to obtain 
a theoretical expression of an average pedestrian flow through the exit. It is found quantitatively 
that the effect of exit door width, a wall, and pedestrian's mood of competition or cooperation 
significantly infiuence the average flow. The results show that there is suitable width of the exit and 
position according to pedestrian's mood. 



I. INTRODUCTION 



Pedestrian dynamics has received growing interest over 
the last decades from physicists since it shows new collec- 
tive behaviors such as dynamicalphase transitions and 
spontaneous symmetry breakingjlf, d, d, 0] ■ Helbing has 
designed the social force model which reproduces typ- 
ical pedestrian behavior such as an arching, lane forma- 
tion, and oscillations of the direction at bottlenecks. It 
is based on a system of coupled differential equations 
which have to be solved by using e.g. molecular dynam- 
ics approach similar to the study of granular matter. An- 
other approach is a discrete modeling using cellular au- 
tomaton, which has been actively studied in recent years 
SEIillil- In this paper, we study the floor field (FF) 
model which is a cellular automaton model, introducing 
two kinds of FFs, i.e.. Static FF (SFF) and Dynamic 
FF (DFF), to move pedestrians from one cell to another. 
The two FFs, which are explained in the Sec|TTl enable us 
to simulate egress processes from complex rooms of arbi- 
trary geometry quite efficiently. Kirchner et al discovered 
that pedestrian's mood of competition or cooperation in- 
crease or decrease the evacuation time [oj. Moreover an 
obstacle in front of the exit will shorten the evacuation 
time in some cases in the simulation of egress processes 
0. Many extended models are proposed up to now to 
make the FF model more realistic. For example, the 
strength of inertia of pedestrians which suppresses quick 
changes of the direction of motion is considered in Q. 
Henein et al has taken into account physical forces be- 
tween pedestrians by adding a dynamic force field to the 
FF model [9]. 

Most of these studies are based on simulations and 
there are few analytical results because of the complex- 
ity of rules of motion and two-dimensionality. In this 
paper we present an analytical result on outflow through 
an exit, which is one of the most important index in eval- 
uating evacuation dynamics. Kirchner et al obtained an 
expression of the average number of evacuated persons 
(TV) from an exit with 1 cell's width as a function of 
time step and the friction parameter ^ by mean-field ap- 
proximation 0]. Here we newly introduce the bottleneck 
parameter [3 which makes pedestrian behaviors around 



the exit more realistic. We have succeeded to calculate 
the average flow {Q) as a function of /3, /it, and the width 
of an exit door w by cluster approximation. As far as 
we know, the analytical expression of the average flow 
through the exit with arbitrary width is derived for the 
first time in this paper. 

This paper is organized as follows. In SecHIlwe briefly 
review the FF model and introduce the new parameter [3 
in Sec lIIII We calculate the average flow {Q) by cluster 
approximation in Sec lIVl and (Q) obtained from sim- 
ulation and the theoretical expression are compared in 
Sec|Vl In Scc lVIl and Sec lVIII we consider how the mood 
of the pedestrians and a wall beside the exit influence 
the average flow. The both effects are explained by the 
contour plots of the average flow in Sec lVIIII Sec|TX] is 
devoted to summary and discussion. 



II. FLOOR FIELD MODEL 

A. Floor Field 

We consider a situation that every pedestrian in a room 
moves to the same exit. The room is divided into cells as 
given in Fig[TJ Man shaped silhouettes represent pedes- 
trians, an alphabet E and alphabets O represent the exit 
cell and obstacle cells, respectively. Each cell contains 
only a single pedestrian at most. Every time step pedes- 
trians choose which cell to move from 5 cells: a cell which 
the pedestrian stands now = (0,0)) and the Neu- 

mann neighboring cells ((i,j) = (0,1), (0,-1), (1,0), 
(—1,0)) (Fig[2]). Two kinds of FFs determine the proba- 
bility of which direction to move. SFF Sij, which is the 
shortest distance to the exit cell, is given by the norm 
as 

Sij — \J — Xexit]'^ + IVij ~ Vexitl'^ : (1) 

where {xij,yij) and {xexit,yexit) are the coordinates of 
the cell (i,j) and the exit cell respectively. However, 
when there is an obstacle on the way to the exit, SFF is 
calculated by making detour of it (FiglSj). Thus, SFF is 
not simply described by ([T]) in the case Q]. Pedestrians 
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move to a cell that has smaller SFF than a cell they oc- 
cupy and hence go to the exit. DFF Dij at cell (i, j) is 
a number of footprints left by the pedestrians. Pedestri- 
ans interact with each other by footprints like ants do by 
their pheromone. The long-ranged interaction by pedes- 
trian's sight is approximated to short-ranged interaction 
around the pedestrian by using DFF. This shortens cal- 
culation time dramatically. When pedestrians move to 
the common exit, it is known in reality that they tend 
to follow each other. This phenomenon can be repro- 
duced by moving pedestrians to a cell which has bigger 
DFF (FigH]). DFF has its own dynamics, namely diffu- 
sion and decay, which leads to broadening, dilution and 
finally vanishing of the footprints [H] . 

Therefore in this model, the transition probability pij 
for a move to a neighbor cell (i,j) is determined by the 
following expression. 



(2) 



Here the values of the FFs Sij and Dij at each cell 
are weighted by two sensitivity parameters and kd with 
the normalization N. There is a minus sign before kg 
since pedestrian move to a cell which SFF decreases, ^ij 
returns for an obstacle or a wall cell and returns 1 for 
other kinds of cells. Note that in our paper a cell occupied 
by a pedestrian is not regarded as an obstacle cell, thus 
it affects the normalization N. 
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FIG. 1: A schematic view of an evacuation simulation by the 
FF model. Pedestrians proceed to the exit by one cell at most 
by one time step. 
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FIG. 2: Target cells for a pedestrian at the next time step. 
The motion is restricted to the Neumann neighborhood in 
this model. 
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FIG. 3: Static floor field constructed by the exit E. The num- 
bers in each cell represent the Euclidean distances from the 
exit cell. 
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FIG. 4: A schematic view of dynamic floor field. At t + 1 two 
pedestrians who could move leave footprints at cells where 
they occupied at t. Remaining pedestrians are likely to move 
cells where footprints are left. 



B. Conflict resolution and friction 

Due to the use of parallel dynamics it happens that 
two or more pedestrians choose the same target cell in 
the update procedure. Such situations are called conflicts 
in this paper. To describe the dynamics of a conflict in a 
quantitative way, friction parameter ^ S [0, 1] was intro- 
duced in Refs.fa, 0|. This parameter describes clogging 
and sticking effects between the pedestrians. In a conflict 
the movement of all involved pedestrians is denied with 
probability ^, i.e., all pedestrians remain at their cell. 
Therefore, the conflict is solved with probability 1 — /i, 
and one of the pedestrians is allowed to move to the de- 
sired cell (FiglS]). The pedestrian which actually moves 
is then chosen randomly with equal probability. In a sit- 
uation with large pedestrians are competitive and do 
not give way to others. Thus they hardly move due to 
the conflict between them. Contrary in a situation with 
small fi they give way and cooperate each other. 



C. Update rules 

The FF model consists of the following 5 steps per unit 
time step, and is repeated until all pedestrians have ex- 
ited or the maximum calculation time steps have passed. 



1. Calculate each pedestrian's transition probability 
by 121) and the values of SFF and DFF. 
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FIG. 5: The way of solving conflicts. In a conflict situation, 
movement of all involved pedestrians remain at their cell with 
probability pt. One of them is randomly allowed to move to 
the desired cell with probability 1 — /i. 



2. Move pedestrians based on the calculated transi- 
tion probability. If there are cells which are possi- 
bly occupied by more than two pedestrians, solve 
conflicts by the means of Sec lII Bl 

3. Diffuse and decay DFF of every cell. 

4. Pedestrians who could move at the step 2 increase 
a value of DFF by 1 at the cell they occupied. 

5. Pedestrians who stand on exit cells are removed 
from the room. 

In the following, we consider only the effect of SFF and 
ignore DFF for simplicity. Thus the step 3 and 4 is not 
considered in this paper. DFF plays a role on mimick- 
ing the long-ranged interaction with pedestrians to short- 
ranged one. Therefore, ignoring DFF is justified when we 
only consider pedestrian behaviors near an exit as done 
here. We also confirm that results in this paper are not 
significantly changed by the introduction of DFF. 



III. INTRODUCTION OF A NEW PARAMETER 
NEAR THE EXIT 



In real situations pedestrian density depends on the 
area in the room. While there are few pedestrians near 
the corner, there are many pedestrians gathering around 
the exit. Therefore, pedestrians often conflict with each 
other around the exit and an arch of pedestrians is likely 
to be formed in front of the exit due to a friction between 
them . Figure [S] shows the number of conflicts in an 
egress process in a competitive situation in 10,000 time 
steps. An exit is set at (a;, y) = (6, 10) in the figure. 
We see 7,842 conflicts at the exit cell, about 1,000 at 
the 5 Moore neighbor cells of the exit, and less than 110 
at other cells. This result says more than 60 percent 
of conflicts occur at the exit cell and the probability of 
conflicting there is about 80 percent. Since pedestrians 
know this fact by experience, they walk fast when they 
are far from the exit, while they walk slowly or give way 
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FIG. 6: The number of conflicts in a competitive situation. 
We simulated in the 11x11 cells' room with the 1 cell's exit 
for 11,000 time steps and accumulated the value of 1,001 to 
11,000 time steps. The cell at (6, 10) is the exit cell. We find 
that more than 60 percent of conflicts occur at the exit cell 
and the probability of conflicting there is about 80 percent. 



to each other around the exit. That is to say, walking 
velocity depends on the area in the room. In the usual 
FF model, however, transition probability, i.e. walking 
velocity is same wherever pedestrians are. To take into 
account this situation, we introduced a new parameter 
(5 £ [0, 1] which we call the bottleneck parameter. The 
transition probability of pedestrians who occupy one of 
the Neumann neighboring cells of the exit cell is described 
as follows: 

= PN^ij exp{-ksSij + kdDij) 

{{^,J)^{0,0)) (3) 
Po,o = il-f3) + PNeM^ksSofl + fcdZ^o.o), 

where N is represented as 

The transition probability of other cells is the same as ([2]). 
Here, if /3 = 0, pij = ^ (0, 0)) andpo^o = 1, which 

means that nobody move to the exit cell. While if /3 = 1, 
the transition probability is the same as ^ , which means 
that pedestrians move fast as they are far from the exit. 
/3 controls the velocity of the pedestrians who are at the 
neighboring cells of the exit. In Ref.[3], the parameter kg 
is used to describe the velocity of the pedestrians. How- 
ever, since small value of kg means lack of a knowledge of 
the exit position, pedestrians sometimes move backward. 
In reality, pedestrians move to the exit along the shortest 
path and slow down near the exit to avoid conflicts with 
others. Therefore, P is not compensated by k^, and we 
expect more realistic pedestrian behavior is seen by the 
parameter /3. When kg is large, the transition probabil- 
ities of pedestrians at neighboring cells of the exit are 
approximated as FiglTl This simplification enables us to 
analyze the pedestrian behavior theoretically. 
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FIG. 7: An approximation of the transition probabilities at 
Neumann neighboring cells of the exit when the width of the 
exit w = 1. When ks is large, po,i — * P, po,o ~* ^ ^ P, and 
po,-i, pi,o, P-1,0 — > for a pedestrian A. Similarly, pi,o — > /3 
for a pedestrian B and p-i,o —> P for a pedestrian C. 



IV. ANALYTICAL EXPRESSION OF THE 
AVERAGE FLOW USING CLUSTER 
APPROXIMATION 



In the /3— introduced FF model. When kg is large, it 
is almost sure that pedestrians move to the exit cells 
with the probability 1 if they are far from the exit and 
with the probability /3 at the Neumann neighboring cells 
of the exit. Therefore, in this section we focus on the 
exit cells and the neighbor cells of them and calculate 
an analytical expression of the average pedestrian flow 
through the exit by cluster approximation. The flow is 
defined as the number of evacuated persons per 1 time 
step thorough an exit. We suppose that a big jam is 
formed around the exit. This enables us to simplify a 
situation that only SFF affect pedestrians' motion. 

First, we calculate the flow when the width of the 
exit w = 1. The transition probability are defined in 
FiglHl We consider two kinds of states of a cell 1 and 
0, which represent that a pedestrian exists at the cell or 
not. Therefore, in the case w ^ \, there are 16 different 
states for these four cells in total. Since we assume a big 
jam exists at the exit, pedestrians enter into three neigh- 
boring cells of the exit with the probability 1. a is the 
probability of getting out from the exit cell, which is set 
as 1 throughout this paper. We define pt(0) as the prob- 
ability that a pedestrian is not at the exit cell at time 
step t and as the probability that a pedestrian is 

at the exit cell at time step t. The master equations are 
described as follows: 
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FIG. 8: Cluster approximation at the exit with one cell. 
a, P2, Ps represent transition probability. 



cells, which is described as follows: 

r = /3i(l - /32)(1 ~ /33) + /32(1 - /?3)(1 - Pi) 
+ /33(l-/3i)(l-/32) 
+ (1 - ^i){Pl|32{l - Ps) + /32/33(l - Pi) 

+ P3Pl{l-P2)+PlP2P3}- (6) 

The first term is the probability of coming a pedestrian 
from the cell A (FiglS]). Similarly, the second and the 
third terms are the probability of coming a pedestrian 
from the cell B and C respectively. The first three terms 
enclosed in parentheses at the forth term represent the 
probability that one of the pedestrians enter into the exit 
cell from two of the three cells (A, B, and C) by resolving 
the confiicts. The last term enclosed in the parentheses 
represents a similar situation, but pedestrians enter into 
the exit cell from all the three neighboring cells. By using 
([5]) and ^ with the normalization condition 



Pti0)+pt{l) = 1 
we obtain the stationary solution 

Poc(l) = 1 



a 



a - a2 - ai - ao - ^{ai + 2ao) 



flo = -P1P2P3 

ai = P1P2 + P2P3 + P3pi 

02 = -(/?! +/32 +/33)- 



(7) 



(8) 



(9) 



Thus the number of pedestrians who can evacuate from 
the exit with one cell's width per a time step, i.e. the 
average pedestrian flow through an exit is described as 
follows: 



{Q{Pi,P2,P3,a,fJ-)) = apoo(l) 



(10) 



The expression of the average flow in the case a = Pi 
P2 = Ps = ^ was obtained in Ref.Q] as 



Here r represents the probability that a pedestrian enter 
into the exit cell from the three Neumann neighboring 



(11) 
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FIG. 9: The three special cases of FigO We assume that 
the arrow with transition probability is interpreted as the 
existence of a wall that blocks pedestrians' motion. 
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which can be recovered by ((5(1, 1, 1, 1, /i)). Thus (jTU]) is 
a generalization of the previous result pip . 

Next we specify pU| by substituting and /3 for (3i, 
/?2, as follows: 



(a) (gi (/?,«)) = {Q{/3, 0,0, a, fi)) = 

(6) {q2W,a,^l)) EE (Q(/3,/?,0,a,/i)) 

a 



1 - 



a + 2/3 - (1 + 
(c) ((?3(/3,a,Ai)) = {Q{P,f3,f3,a,fi)) 

a 



a 



1 



a + 3/3 - 3(1 + Ai)/32 + (1 + 2/i)/33 



These expressions describe average flows through an exit 
with the configuration described in Fig[5] respectively. 

Finally we calculate the average flow of pedestrians 
through an exit with arbitrary w S N width. The re- 
lation between the width of an exit and the outflow is 
an important study which has been investigated experi- 
mentally so far [10, . When pedestrians move to the 
exit cell along the shortest path, social morals may sup- 
press pedestrians breaking into the line. Thus, they do 
not gather around the exit in disorder, but tend to form 
lines in front of the exit. Moreover they do not easily 
change lanes in a crowd situation. There are also the ex- 
perimental results that the pedestrian outflow increases 
linearly as the width of an exit increases [ll[. Therefore, 
we can represent the average flow through the exit with 
w cell's width by linear sum of {qi), ((72), and (qs). Here 
we consider two types of exits: an exit at the center of 
the wall (Ce-exit) and an exit at the corner of the room 
(Co-exit). Ce-exit {w > 2) is divided into (gi)-exits, and 
(q'2)-exits, and an average flow through an exit with w 
cell's width {Q center, w) IS described as: 



(0 



center,w 



(93) 



(u;-l) 



2{q2) + {w~2){qi) {w>2). 



(13) 



In a similar way, Co-exit is divided into (gi)-exits, and 
(q2)-exit, and an average flow through an exit with w 
cell's width {Qcorner,w) IS described as follows: 



FIG. 10; Dividing an exit with 3 cells into three exits with 1 
cell. The outflow is also represents the sum of each flow. 



The examples of dividing Ce-exit and Co-exit {w = 3) 
are shown in FigfTUl We also define the average flow per 
1 cell as: 



(12) 



{Qcenter) {Qcenter^w) /'^ 
{qcorner) = {Q corner,w )/w. 



(15) 



corner.w 



) = (<Z2) + (w-i)(gi) i'w>i). 



(14) 



V. COMPARING THE AVERAGE FLOWS OF 
ANALYSIS AND SIMULATION 

In this section we compare the analytical and computa- 
tional results of the average flow (q) , which is a function 
of p. The parameters of FFs are set at: kg = 10 and 
kd — 0.0 so that pedestrians move to the exit along the 
shortest path. The size of the room used in simulation is 
11x11 cells. Ce-exit room has the exit cells at the center 
of the boundary cells of one side of the room, and those 
of the other sides are all entrance cells where pedestrians 
come with the probability 1. Similarly, Co-exit room has 
the exit cells at the corner of the room, and boundary 
cells of two sides which do not include the exit cells are 
all entrance cells. The examples of the 5x5 rooms are 
shown in Fig[TlJ We simulated 11,000 time steps with 
the initial condition that pedestrians occupy all cells ex- 
cept exit and obstacle cells. Then the average flow of 
10,000 time steps from 1,001 to 11,000 is used to depict 
FigdU It shows average pedestrian flows at the exit as 
a function of /3 for various /i values. We see that the 
simulations agree with the analytical results very well. 
The errors become large for /i = 0.9 since pedestrians 
conflict with each other and then they cannot move not 
only at the exit cells but also at the other cells. Sur- 
prisingly, for /I = 0.9, we clearly find that a maximum 
fiow is attained at a value of /3 in both simulations and 
analytical results. We call it as the optimal /3, notated as 
Popt here after. For /i — > 1 the number of unsolved con- 
fiicts increases as /3 grows. As a result, pedestrians stick 
and the average pedestrian flow decreases. We also find 
that the differences of the flows for different /i is getting 
smaller as w increases by comparing the lengths of the 
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FIG. 12: Average flow (g) as a function of /3 for different [i, 
exit width w, and exit location. (a)Ce-exit, w = \. (b)Co- 
exit, w = \. (c)Ce-exit, w = "i. (d)Co-exit, w = 3. For 
/i = 0.9 we clearly see a maximum of flow at an optimal 
/3. The lengths of the arrows A, B, C and D represent the 
differences of the flows for different ^. 



arrows in FigfT2l i.e., the arrow C is shorter than A and 
D is shorter than B. For the same w the differences of 
the flows for different /x are smaller at the Co-exit than 
the Ce-exit since the arrow B is shorter than A and D is 
shorter than C. 



VI. COMPETITIVE, COOPERATIVE 
BEHAVIOR AND THE WIDTH OF AN EXIT 

Kirchner et al made a research on how the pedestrians' 
mood influences on the evacuation time Q. Both the 
experimental and the computational results show that 
competition is beneficial if the exit width exceeds a cer- 
tain width, and harmful if the exit width is lower than it. 
We explain this phenomenon by the analytical solutions 
(Uni) and (HI]). In Ref.j^, competition is described as a in- 
creased assertiveness (large k^) and a strong hindrance in 
conflict situations (large ji). Cooperation is represented 
by small ks and vanishing /x. In our new Model we de- 
scribe the assertiveness by p. We use /z as a parameter 
of hindrance in conflict situations and its values are same 
as in Ref.Q. The parameters are set at f3 = 1.0, fi = 0.6 



FIG. 13: Average flow for various exit door width w. (a) Ce- 
exit. (b)Co-exit. We observe the crossing of the two curves 
in both (a) and (b). 



for the competitive situation and f3 = 0.4, /i — 0.0 for 
the cooperative situation. 

Figure [T3] shows the average flows for variable door 
width w. The results of analysis agree with those of sim- 
ulation very well. The simulation condition is the same 
as we described in SeclVj however, we used 12x12 rooms 
to set up the exit at the center of the boundary cells of 
the room if the width of the exit is even number. The 
size of the room does not influence the average flow since 
most of the cells are occupied by pedestrians. Clearly we 
can observe the crossing of the two curves at a critical 
door width ~ 3 in Fig |13( aV Our result is well cor- 
respond to the Ref. ^]'s result, which is Wc ~ 2.5. This 
means that we should cooperate with each other to in- 
crease the average pedestrian flow when the width of the 
exit is narrow. On the contrary, when the width of the 
exit is wide, we do not have to give way to other pedes- 
trians and should go through the exit aggressively. When 
the exit door is at the corner of the room, the crossing is 
observed at Wc ~ 1.5 in Fig fT5f bV Therefore, Co-exit is 
more suitable for competitive situation than Ce-exit. 

The Japanese building standards law [l6| gives an 
average pedestrian flow 1.5[persons/(m-s)] if an exit is 
directly connected to the ground. We find that this 
value significantly changes by the pedestrians' moods, 
i.e., competitive and cooperative. From Fig |13l we 
obtain the values of the average flow through Ce- 
exit i.e., 1.5[persons/(m-s)] in the competitive situa- 
tion and 2.0[persons/(m-s)] in the cooperative situation. 
The values are calculated by defining the cell size as 
50[cm] x50[cm] and using a pedestrian velocity 1.3[m/s] 
which is according to the Japanese building standards 
law. 



VII. COMPETITIVE, COOPERATIVE 
BEHAVIOR AND THE EFFECT OF A WALL 

Here we compare the average fiows of Ce-exit and Co- 
exit, and discuss how the wall has an effect on them. The 
difference of {Qcenter,w) and {Qcorner,w) is calculated as 



FIG. 14: The curves of (3^ on /?— plane. (a)io = 1. (b)w > 2. 
The flow of Co-exit is larger than the flow of Ce-exit in the 
upper right region in the flgures. In the lower left region, the 
Ce-exit flow is larger. 



follows: 

{Q corner, w ) - {Qcenter,w) 

^ f-/3(/3-TT2^)(/?-lM (^ = 1) (16) 
\PiP-TTJl)B {w>2), 

where A and B are positive in the entire domain of /? and 
II that is described as 



A 



a2(l + 2^) 



{a + 2/3- (l + ^)/32} 



B = 



{a + 3f3- 3(1 + Ai)/32 + (1 + 2/i)/33} 
a2(l-f 



(17) 



{a + p){a + 2p-{l + ii)P^}' 



We obtain Pc which is the value of /? that {Q center, w) 
equals {Qcorner.w) as foUows: 



/3c = 




{w = l) 
{w > 2), 



(18) 



The curves of p8|) are drawn in Fig[T4l They divide the 
(3 — fj, plane into two regions. In the lower left region, the 
Ce-exit flow is larger and in the upper right region the 
Co-exit flow is larger. We also plot the competitive and 
cooperative situation used in Sec lVIl The figures show 
that the Co-exit flow is larger in the competitive situation 
since the wall prevents pedestrians rushing to the exit 
at the same time, but the Ce-exit flow is larger in the 
cooperative situation. From this result, we can say that 
an exit should be at the center of a wall when pedestrians 
are in the cooperative mood, and should be at the corner 
of the room when people are in the competitive mood. 
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FIG. 15: Contours of average flows for various width of the 
exit and its position. (a)Ce-exit, «; = 1. (b)Co-exit, w = 1. 
(c)Ce-exit, w = 3. (d)Co-exit, w = 3. represents compet- 
itive situation and represents cooperative situation given in 
Sec lVIl Black bold curves represent /3opt, which move to the 
right in order (a), (b), (c), and (d). In Fig. (a) and (b) flow 
in the cooperative situation is larger than that of the com- 
petitive situation. In Fig.(c) the average flows are almost the 
same in both the competitive and the cooperative situation. 
In Fig.(d) flow in the competitive situation is larger. These 
figures explain that flow in the competitive situation is getting 
larger than that in the cooperative situation by increasing the 
width of the exit and the effect of wall. 



we have set at a = 1 in this paper, the average flow is 
determined by /3 and /i. FiglTS] are contour plots of the 
average flow, in terms of (3 and fj,. FigHni (a), (b), (c), 
and (d) correspond to Fig[T2](a), (b), (c), and (d) respec- 
tively. The flow is large in the white region and small in 
the black region. The values of the flow are normalized 
in each figure for drawing the gray scaled figures. The 
thick curves in the figures are the curves of Popt, which 
gives a maximum average flow in the case of a constant 
fJ" f^opt curves divide the plane into two regions. In the 
upper right regions pedestrians should slow down further 
to avoid a conflict and in the lower left regions they have 
to speed up to increase the average flow around the exit. 
We calculate /3opti, /3opt2, and (3opt3 from (qi), ((72}, and 
{qs) given in (fT2|l . respectively as follows: 



VIII. A CHANGE OF CONTOUR PLOTS OF 
THE AVERAGE FLOW 

The average pedestrian flow through an exit is decided 
by three parameters: a, (3, and according to pop . Since 



opt! 



1, 



1 



'opt 2 



optZ 



l + Ai' 
1 

1 + 2^' 



(19) 
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FIG. 16; The curve of Popti, (3opt2, and Popts- We see that 
the curve of popta is at the most left, that of f3opt2 is in the 
middle, and that of /3opti is at the most right. 



The Popt of {Q center, w) and {Q corner, w) are straight- 
forward, but not expressed in a simple form, so we omit 
them in this paper. We also plot the competitive and 
cooperative situations used in Sec lVIl in the figures. 

We see that the f3opt curves move to the right in 
the order: (a), (b), (c), and (d). This is explained by 
([ni)C([Tll)Cand First, we compare the (iopt curve's 

position of (iopti, Popt2, and f3opt3 in the /3 - ^ plane by 
P^ . The curve of Popts is at the most left, that of Popt2 
is in the middle, and that of Popti is at the most right 
fFig fTH)) . Next, the expressions of the average flow corre- 
sponding to (a), (b), (c), and (d) are described as follows: 



Q{a) = 


{Q center ("^ 


-1)) 


= (93) 




Qib)^ 


i^Q corner (^U) 


-1)) 


= (92) 






{Q center ('^ 


-3)) 


= 2(92) 


^(91) 


Q{d) = 


i^Q corner (^U) 


-3)) 


= (92) 4 


-2(gi) 



(20) 



Now we see clearly why the f3opt curve of Q{b) is at 
more right than that of Q{a). Q{c) includes ((72) + {qi) 
more than Q{b), and Q{d) includes 2(gi) more than (5(6). 
Therefore, j3opt curve moves to the right in order (a), (b), 
(c), and (d). 

In the lower left regions of the figures the average flow 
increases as /3 increases, but in the upper right regions 
it decreases as (3 increases for fixed ^. Thus the fiow in- 
creasing region expands as Popt curves move to the right. 
This makes an exit more suitable to the competitive sit- 
uation than the cooperative situation. We can interpret 
that increasing of the width of a exit and the effect of 
a wall make the average flow larger in competitive situ- 
ation than cooperative situation, since (3opt curves move 
to the right by both increasing of the width of a exit 
((a)(c)C(b)(d)) and the effect of a wall ((a)(b)C(c)(d)). 



IX. CONCLUSION 

We have introduced it to the FF model that the ef- 
fect of slowing down of pedestrians around an exit, and 
obtained the analytical expression of the average flow 
through an exit with arbitrary w cells by employing clus- 
ter approximation. It turns out that the theoretical re- 
sults agree quite well with the simulations. The effects 
of pedestrians' mood, a width of an exit, and wall ef- 
fect are also studied. The critical exit door width, which 
was first obtained experimentally and was reproduced by 
simulations in Ref. 6] , is also analytically obtained in this 
paper. We find that an exit should be at the center of a 
wall in the cooperative situation whereas it should be at 
the corner of the room in the competitive situation for 
smooth evacuation. The theoretical results also tell us 
that the unsolved conflicts between pedestrians around 
the exit are the main cause of decrease in the average 
pedestrian fiow. Therefore, we should consider how to 
decrease confiicts at a bottleneck to get large pedestrian 
outfiow. 

It is important to study pedestrian's behavior quanti- 
tatively by theoretical analysis, since its dynamics are 
mainly studied by simulations so far. The Japanese 
building standards law gives the average pedestrian flow 
through an exit as a constant value 1.5[persons/m-s] [ig} . 
Our expression of the average pedestrian fiow is more pre- 
cise and realistic, thus our results can be applied to the 
design of buildings so that pedestrians evacuate safely 
and quickly. For example, many present concert halls 
have an exit at the center of the wall, however, accord- 
ing to our study we can shorten an evacuation time by 
setting up an exit at the corner of the hall when people 
rush into the exit in competitive mood. 

In this paper we consider the average fiow through the 
exit with more than 1 cell as the linear sum of the flow 
through an exit with 1 cell. In the calm situation, social 
morals suppress pedestrians cutting into lines, however, 
in the panic situation, they break into lines to save their 
lives. Introducing such interactions between neighboring 
cells of an exit in detail is the future works. 
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